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^ I Abstract 

^-j' In this paper, results are presented on studies which have been performed to investi- 

gate the impact of the cut-off elevation angle (CEA) and elevation-dependent weighting 
O ' (EDW) on the EOP estimates and baseline length repeatability. For this test, CONT05 

observations were processed with different CEA and EDW, keeping all other options 
the same as used during the routine processing. Uncertainties and biases, as well as 
correlations between estimated parameters have been investigated. It has been shown 
that small CEA, up to about 8-10 degrees does not have large impact on the results, 

y—i . and applying EDW allows us to get better result (smaller errors). However, this result 

has been proven with standard geodetic VLBI observations, where rather few observa- 

^ | ' tions were made at low elevations. Perhaps, special R&D sessions with more uniform 

distribution of observations over elevation may be useful for more detailed study on 
the subject. 
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1 Introduction 



It is well known that precision and accuracy of astronomical observations, both optical and 
radio, made through the Earth's atmosphere depend on the elevation at which the object 
is observed. These errors grow with decreasing of the elevation due to larger air mass 
and difficulties in modelling of refraction effects at low elevation. From this point of view 
observations should be made in the near-zenith zone when possible. 

On the other hand, inclusion in processing of observations made at low elevations is 
important when definite groups of highly correlated parameters, for instance station coor- 
dinates and zenith troposphere delays, are estimated simultaneously. In such a case using 
observations made at in a widest range of elevation allows us to mitigate the correlations 
between unknowns and improve the solution. 

To meet these mutually exclusive requirements, proper elevation-dependent weighting 
(EDW) of observations is used. In a special case of step-like weighting function, i.e. rejection 
of the observations made at the elevation less than the given limit, such a limit usually is 
called cut-off elevation angle (CEA). 

It was shown in many studies that elevation-dependent weighting may have a significant 
impact on the results of processing of the space geodesy observations. In particular, several 
studies of this effect was made by the Goddard and Vienna VLBI analysis groups in the 
framework of the IVS VLBI2010 Committee activit}0. They investigated an influence of 
CEA and EDW on geodetic results such as Earth orientation parameters (EOP), baseline 
length repeatability, troposphere parameters, and station heights. Those results were based 
on simulation. Gipson in [TJ used another approach to EDW. He applied elevation dependent 
additive noise to the measurement error instead of using a weighting factor as it is usually 
being made. He tested his method with the actual CONT05 observations. Results of both 
mentioned and other results are sometimes contradictory. This gave an impulse to the 
present work where some results are presented of investigation of the impact of the CEA 
and EDW on the baseline length repeatability and EOP estimates. 

2 Test description 

For this test, CONT05A observations were processed making use of OCCAM software with 
different EDW functions including CEA, keeping all other options as follows: 

• Kalman filter mode (KF), 

• random walk model for clocks, PSD=1.5 ps 2 /s, 

• random walk model for ZTD, PSD=0.25 ps 2 /s, 

• one NS and EW troposphere gradient estimate for the session. 

For the continuous EDW mode (continuous weighting function), the measurement error 
is multiplied by a factor 

W e = (sine /sine) p , (1) 
1 href =" http : / /ivscc . gsfc . nasa. gov /pipermail/ ivs- v2c / 
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Figure 1: Actual values of the hydrostatic (HMF), wet (WMF) mapping functions for the 
C0509 (05SEP20XA) session (one point correspond to one observation). Solid line corre- 
sponds to function 1/ sine. Data for the full elevation range (left) and zoomed data for low 
elevations (right) are shown. 



where en and p are EDW parameters, e is the source elevation. Such a weighting function 
provides a smooth stepless change in weight for any en. One can see that a case of en = 
90°, p — 1 gives merely W e = 1/ sine, which is close to actual wet mapping function used 
in the last works by MacMillan and Gipson (private communications). Figure [T] shows 
actual hydrostatic and wet mapping functions along with approximation 1 / sin e for a typical 
CONT05 session. One can see that these functions are close enough, and either of them can 
be used for the EDW without significant impact on the result. 

EDW mode with e = 10, p = 2 was implemented in the OCCAM/GROSS software [2] 
for routine data processing. It will be referred hereafter as "normal mode". 

In a case of CEA we have 

w = | 1, if e >= e 
1 10 3 otherwise 

The latter line corresponds to the KF realization used in OCCAM. 

For VLBI delay, measurement error coming from correlator is multiplied by two W e 
values computed for both the stations. In our test, en = 3(2)25° were used for CEA test, 
and en = 10, 25, 45, 90°, p — 1, 2 were used for continuous EDW mode. 



3 Test Results 

3.1 Baseline length repeatability 

Test results obtained with different CEA are shown in Figure [2j The case of en = 3° includes 
all the observations without weighting, since no CONT05 observations were made at the 
elevation less than 4°. 

Table [T] shows EDW test results. Different EDW modes are denoted as w_e_p, where e 
and p are e and p in Eq.([T]). Test results are given for quadratic approximation in percent 
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Figure 2: Dependence of the baseline length repeatability on the cut-off elevation angle 
(quadratic regression): all tested modes (left) and data for low elevations (right). 
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Table 1: Comparison of the baseline repeatability obtained with different EDW models. See 
explanation in text. 



EDW mode 


Baseline len 


gth, 10 3 km 




3 


6 


9 


12 


w_fO_f 


98.0 


97.2 


97.3 


97.6 


w_f0_2 


95.3 


95.0 


95.5 


96.1 


w_25_f 


102.6 


92.8 


87.2 


83.8 


w_25_2 


95.7 


91.8 


91.4 


91.9 


w_45_f 


101.7 


90.0 


84.6 


81.9 


w_45_2 


104.0 


97.2 


101.8 


107.6 


w_90_l 


101.3 


90.5 


85.4 


82.8 


w_90_2 


109.6 


127.3 


161.1 


191.1 



Table 2: EOP statistics for different EDW models. See explanation in text. 



Statistics 


All 


EDW mode 






10_1 


10_2 


25_1 


25_2 


45_1 


45_2 


90_1 


90_2 


Xp uncertainty 


26 


26 


25 


26 


28 


27 


37 


30 


51 


Yp uncertainty 


25 


24 


24 


24 


26 


25 


32 


26 


43 


UT1 uncertainty 


1.1 


1.1 


1.1 


1.1 


1.1 


1.0 


1.3 


1.1 


1.6 


Xc uncertainty 


21 


21 


20 


20 


19 


20 


22 


22 


27 


Yc uncertainty 


19 


18 


18 


18 


18 


18 


21 


20 


27 


Xp bias w.r.t. IGS 


-79 


-75 


-79 


-77 


-73 


-80 


-86 


-81 


-123 


Yp bias w.r.t. IGS 


+135 


+135 


+128 


+132 


+129 


+134 


+127 


+129 


+92 


Xp wrms w.r.t. IGS 


77 


72 


76 


69 


65 


68 


69 


64 


86 


Yp wrms w.r.t. IGS 


73 


74 


72 


69 


61 


68 


63 


71 


83 


mean 


75 


73 


74 


70 


63 


68 


66 


68 


84 



with respect to the case of CEA with e$ = 3°. One can see that several EDW modes show 
about the same improvement in the baseline length repeatability. 

3.2 EOP 

Figure [3] shows EOP statistics for different CEA. Notation used is the following: Xp, Yp - 
terrestrial pole coordinates, Xc, Yc - celestial pole coordinates. Weighted Allan deviation 
is computed as described in [3]. All the results related to Xp, Yp, Xc and Yc are given in 
^as, the results related to UT1 are given in /is. Comparison with IGS time series is shown 
in Figure HI Xp and Yp wrms with respect to IGS EOP series are computed after removing 
the bias. 

Table [2] shows the main EDW test results. Notation of EDW modes is the same as in 
Table [1] with the 'w_' prefix omitted. 'All' column correspond to inclusion of all observations 
without weighting. One can see again that several EDW modes show about the same EOP 
precision and accuracy. 
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Figure 3: Statistics of EOP obtained from CONT05. 
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Figure 4: Comparison of EOP obtained from CONT05 with the IGS EOP series. 
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4 Conclusions 



The preliminary conclusions from this test are the following. 

• The baseline length repeatability steadily grows with the CEA increasing, remaining 
practically the same in the cut-off angle range from 3° (i.e. no cut-off for the CONT05) 
to 9°. 

• The best result is obtained when the EDW elevation-depending weighting is applied to 
the low-elevation observations. However, the test results are not always unambiguous. 
Further adjustment of the weighting method may be fruitful. 

• The Xp, Yp and UT1 uncertainties grow with the increasing cut-off angle after about 
10°. Most probably, this reflects the fact that only about 6% of the total number 
of CONT05 observations were made at the elevations below 10°. The Xc and Yc 
uncertainties and scatter depend on the CEA much less. 

• Xp bias w.r.t. IGS slightly depends on the CEA, except the maximum tested CEA 
values, evidently unrealistic. In contrast, Yp bias substantially changes with increasing 
CEA. Most probably, this can be explained by the CONT05 network orientation, for 
which the longitude of the central meridian Ao = 265° just corresponds to the Y 
direction of the terrestrial coordinate system. 

• Some statistics such as the uncertainty and the scatter of the Xc and Yc, as well as the 
WRMS of Xp and Yp w.r.t. IGS have the minimum at the CEA around 15°, which is 
interesting and deserves a supplement investigation. 

• As one can expect, the correlations between EOP comprising Xp and Yp grow with 
increasing CEA, but remain small due to good CONT05 network geometry. The same 
can be expected for the IVS2010 network. The correlation between Xc and Yc re- 
main practically the same for all tested CEA, except the maximum tested CEA value, 
evidently unrealistic. 

Finally, we can conclude that inclusion of the low-elevation observations, properly weighted, 
improves the baseline length repeatability and EOP results. On the contrary, filtering the 
observations using the cut-off elevation method may lead to degradation of geodetic results. 
However, this should be mentioned that the conclusions drawn from the result obtained in 
this paper has been proven with standard geodetic VLBI observations, where rather few ob- 
servations were made at low elevations, as mentioned above. Perhaps, special R&D sessions 
with more uniform distribution of observations over the sky, including observations at very 
low elevations, may be useful for more detailed study on the impact and optimal processing 
of the low-elevation observations on geodetic parameters obtained from VLBI observations. 

It ought be mentioned that all the EDW modes considered in this paper in fact modify 
only diagonal elements of the corresponding covariance matrix. According to Gipson's work 
P] best result can be achieved in a case of account also for correlations between observations. 
It seems to be interesting to investigate how this approach will work in KF estimator. 
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